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Abstract 

We obtain an expression for the active gravitational mass of a rela- 
tivistic heat conducting fluid, just after its departure from hydrostatic 
equilibrium, on a time scale of the order of relaxation time. 

It is shown that an increase of a characteristic parameter leads 
to larger (smaller) values of active gravitational mass of collapsing 
(expanding) spheres, enhacing thereby the instability of the system. 

1 Introduction 

In a recent series of works Jl], ^ [| the behaviour of dissipative systems at 
the very moment when they depart from hydrostatic equilibrium, has been 
studied. 

It appears that a parameter formed by a specific combination of thermal 
relaxation time, temperature, proper energy density and pressure, may crit- 
ically affect the evolution of the object. 
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More specifically, it has been shown that in the equation of motion of any 
fluid element, the inertial mass term is multiplied by a factor, vanishing for 
a given value of that parameter (critical point) and changing of sign beyond 
that value. 

Although the above mentioned parameter is constrained by causality require- 
ments, it appears that in some cases these requirements do not prevent the 
system from reaching the critical point ||. Furthermore it might not be 
reasonable to apply, close to the critical point, restrictions obtained from a 
linear perturbative scheme (as is the case for causality conditions) 0. 
In order to delve more deeply into the physical nature of the critical point, 
we shall obtain here an expression for the active gravitational mass (Tolman 
mass which explicitly contains the parameter mentioned above. 
It will be seen that this expression yields larger (smaller) values for the active 
gravitational mass of the inner core of a collapsing (expanding) sphere, as we 
approach the critical point, this tendency persists beyond the critical point 
as the systems moves away from it. 

This result provides some hints about the way in which the evolution of the 
system is affected by the aforesaid parameter. 
The paper is organized as follows. 

In the next section the field equations, the conventions and other useful 
formulae are introduced. In section 3 we briefly present the equation for the 
heat conduction. The departure from hydrostatic equilibrium is considered 
in section 4. In section 5 we derive an expression for the Tolman mass 
and evaluate it at the very moment the system departs from hydrostatic 
equilibrium . Finally a discussion of this expression is presented in the last 
section. 

2 Field Equations and Conventions. 

We consider spherically symmetric distributions of collapsing fluid, which for 
sake of completeness we assume to be anisotropic, undergoing dissipation in 
the form of heat flow, bounded by a spherical surface S. 
The line element is given in Schwarzschild-like coordinates by 

ds 2 = e u dt 2 - e x dr 2 - r 2 (d6 2 + sin'Odcf 2 ) (1) 
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where v{t,r) and X(t,r) are functions of their arguments. We number the 

coordinates: x° = t; x 1 = r; x 2 = 8; x 3 = 0. 

The metric (]!]) has to satisfy Einstein field equations 

Gl = -SttT; (2) 

which in our case read f5|: 

_ 8l r» = -I + e-(i-^) (3) 



- 8ttT 2 2 = -8ttT| = - — (2A + A(A - u 

+ ^(2," + ,' 2 -AV + 2^) (5) 

- 8ttT 01 = -- (6) 
r 

where dots and primes stand for partial differentiation with respect to t and 
r respectively. 

In order to give physical significance to the components we apply the 
Bondi approach ||. 

Thus, following Bondi, let us introduce purely locally Minkowski coordinates 

(r,x,y,z) 



dr = e v l 2 dt dx = e x ^ 2 dr dy = rd9 dz = rsin 

Then, denoting the Minkowski components of the energy tensor by a bar, we 
have 



01 
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Next, we suppose that when viewed by an observer moving relative to these 
coordinates with proper velocity u in the radial direction, the physical con- 
tent of space consists of an anisotropic fluid of energy density p, radial pres- 
sure P r , tangential pressure P± and radial heat flux q. Thus, when viewed 
by this moving observer the covariant tensor in Minkowski coordinates is 



P 

-Q 






-Q 

Pr 










P± o 

P, 



/ 



Then a Lorentz transformation readily shows that 



rp{] 

J- n 



-0 _ p + P r UJ 2 

J ~~ - n "T 



T\ = T\ = 



P r + puj" 1 



2Qcue 



A/2 



(i -u 2 yi 2 



2Que 



A/2 



UJ 



!)V2 



± 



(7) 

(8) 
(9) 



loi-e 1 01 - (1 - uj 2 ) 1 / 2 ^ ' ^ ' 

with 

Note that the coordinate velocity in the (t,r, 9, (ft) system, dr/dt, is related 
to uj by 

w = ^: e (A-)/2 (12) 

dt 

At the outside of the fluid distribution, the spacetime is that of Vaidya, given 
by 

ds z = (i_ du 2 + IdudlZ - K 2 (d6 2 + sin 2 9d<t> 2 ) (13) 
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where u is a time-like coordinate such that u = constant is (asymptotically) 
a null cone open to the future and 7Z is a null coordinate (gnu = 0). It 
should be remarked, however, that strictly speaking, the radiation can be 
considered in radial free streaming only at radial infinity. 
The two coordinate systems (t, r, 9, 0) and (u, 7Z, 9, <fi) are related at the 
boundary surface and outside it by 

u = t-r - 2Mln(— - l) (14) 

71 = r (15) 

In order to match smoothly the two metrics above on the boundary surface 
r — r T,{t), we have to require the continuity of the first fundamental form 
across that surface. As result of this matching we obtain 



[Pr] 



Q e A / 2 fl 
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1/2 



(16) 



expressing the discontinuity of the radial pressure in the presence of heat 
flow, which is a well known result |3|. 

Next, it will be useful to calculate the radial component of the conservation 
law 

= (17) 
After tedious but simple calculations we get 



(-*< = ^ (r; - if) + W (r; - is) + tl { X + x - - x ± ) (is) 



which in the static case becomes 

representing the generalization of the Tolman-Oppenheimer-Volkof equation 
for anisotropic fluids 0. 
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3 Heat Conduction Equation. 



In the study of star interiors it is usually assumed that the energy flux of 
radiation (and thermal conduction) is proportional to the gradient of tem- 
perature (Maxwell- Fourier law or Eckart-Landau in general relativity). 
However it is well known that the Maxwell-Fourier law for the radiation flux 
leads to a parabolic equation (diffusion equation) which predicts propagation 
of perturbation with infinite speed (see |§-[1(| and references therein). This 
simple fact is at the origin of the pathologies [TTI] found in the approaches of 
Eckart [12] and Landau Jl3| for relativistic dissipative processes. 
To overcome such difficulties, different relativistic theories with non- vanishing 
relaxation times have been proposed in the past |L1|-||17||. The important 



point is that all these theories provide a heat transport equation which is 
not of Maxwell- Fourier type but of Cattaneo type [113], leading thereby to a 



hyperbolic equation for the propagation of thermal perturbation. 
Accordingly we shall describe the heat transport by means of a relativistic 
Israel-Stewart equation [|l(J. Although a complete treatment of dissipative 



processes requires the inclusion of viscous stresses as well as the coupling 
between these and the heat flow, we shall assume here for simplicity vanishing 
viscosity. Thus we have 



t ^d7 + qa = KpaP (T> " Tap) ~ TU ° q ^ ~ \ kT2 {^ u ^ - qa (20) 



;/3 



with 



e-^ 2 uoe- x ' 2 



q* = Q (ue {x - v)/ \ 1, 0, 0) (22) 

where k, r, T, q 13 and aP denote thermal conductivity, thermal relaxation 
time, temperature, the heat flow vector and the components of the four 
acceleration, respectively. Also, P a @ is the projector onto the hypersurface 
orthogonal to the four velocity u a . 
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4 Thermal Conduction and Departure from 
Hydrostatic Equilibrium. 

Let us now consider a spherically symmetric fluid distribution which initially 
may be in either hydrostatic and thermal equilibrium (i.e. to = Q = 0), or 
slowly evolving and dissipating energy through a radial heat flow vector. 
Before proceeding further with the treatment of our problem, let us clearly 
specify the meaning of "slowly evolving". That means that our sphere 
changes on a time scale which is very large as compared to the typical time 
in which it reacts on a slight perturbation of hydrostatic equilibrium. This 
typical time is called hydrostatic time scale. Thus a slowly evolving system 
is always in hydrostatic equilibrium (very close to), and its evolution may be 
regarded as a sequence of static models linked by (|6]). This assumption is 
very sensible, since the hydrostatic time scale is usually very small. 
In fact, it is of the order of 27 minutes for the sun, 4.5 seconds for a white 
dwarf and 10 -4 seconds for a neutron star of one solar mass and 10 Km 



radius [19 



In terms of uo and metric functions, slow evolution means that the radial 
velocity u measured by the Minkowski observer, as well as time derivatives 
are so small that their products and second order time derivatives may be 
neglected (an invariant characterization of slow evolution may be found in 

@)- 
Thus IE! 



i)^\^Xu^X 2 ^u 2 ^uj 2 ^dj = (23) 

As it follows from (El) and (10), Q is of the order 0(u). Therefore in the 
slowly evolving regime, relaxation terms may be neglected and ( p0|) becomes 
the usual Landau-Eckart transport equation [^l[ . 

Then, using (^) and ([T^) we obtain flT9|), which as mentioned before is the 
equation of hydrostatic equilibrium for an anisotropic fluid. This is in agree- 
ment with what was mentioned above, in the sense that a slowly evolving 
system is in hydrostatic equilibrium. 

Let us now return to our problem. Before perturbation, the two possible 
initial states of our system are characterized by: 

1. Static 

Co = Q = oj = Q = (24) 
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2. Slowly evolving 

uj = Q = (25) 
Q « 0(w) ^ (sma/Z) (26) 



where the meaning of "small" is given by ( |23|) . 

Let us now assume that our system is submitted to perturbations which force 
it to depart from hydrostatic equilibrium but keeping the spherical symmetry. 
We shall study the perturbed system on a time scale which is small as com- 
pared to the thermal adjustment time. 

Then, immediately after perturbation ("immediately" understood in the 
sense above), we have for the first initial condition (static) 

u = Q = (27) 

Co « Q jL {small) (28) 
whereas for the second initial condition (slowly evolving) 

Q « 0(u) ^ (small) (29) 

Q w Co ^ (small) (30) 

As it was shown in M, it follows from ( TTH ) and (|7|)-(|I0D that after perturba- 
tion, we have for both initial conditions (see |1] for details) 

- e (-A)/ 2jR = (p + ^ + g e A/2 (31) 

where P denotes the left-hand side of the TOV equation, i.e. 



^ _ dP r 47rrP r 2 i P r m i AnrpP, 



dr l — 2m/r r 2 (l — 2m/r) 1 — 2m/r 
pm 2 (Pj_ - P r ) 



r 2 (1 — 2m/r) r 

P'r + j(P + P r) - 2 -( P t- P r) (32) 
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The physical meaning of R is clearly inferred from (0). It represents the 
total force (gravitational + pressure gradient + anisotropic term) acting on 
a given fluid element. Obviously, R > 0/R < means that the total force is 
directed inward /outward of the sphere. 



Let us now turn back to thermal conduction equation fl20|) . Evaluating it im- 
mediately after perturbation, we obtain for both initial configurations (static 
and slowly evolving) (see for details) 

T Qe x/2 = -kTCj (33) 
Finally, combining (|3ll) and (0), one obtains 

1 

(" + f » (1-^) 
or, defining the parameter a by 

kT 

r(p + P r ) V ; 

- e W 2 R= (p + P r )u(l-a) (36) 

This last expression has the obvious "Newtonian" form 

Force = mass x acceleration 
since, as it is well known, (p + P r ) represents the inertial mass density and by 
"acceleration" we mean the time derivative of u and not (a^a^) 1 / 2 . If a < 1, 
then an outward /inward acceleration {to > 0/u < 0) is associated with an 
outwardly / inwardly (R < 0/R > 0) directed total force (as one expects!). 
However, if a = 1, we obtain that Cj ^ even though R — 0. Still 
worse, if a > 1, then an outward /inward acceleration is associated with 
an inwardly / outwardly directed total force!. 

As mentioned before, the critical point may be restricted by causality condi- 
tions, particularly in the pure bulk or shear viscosity case however this 
is not so in the general case [[|. 

Independently of this fact, it is clear from PB]), that the "effective" inertial 
mass term decreases as a increases. In the next section we shall obtain an 
expression for the active gravitational mass explicitly containing a. 
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5 The Tolman mass 



The Tolman mass for a spherically symmetric distribution of matter is given 
by (eq.(24) in 1) 



rriT 



+ 



4tt rV^ 2 (T ° - Tl - 2T|) dr 
d ( dL 



2 Jo 



r 2 e {u+X)/2 



dt \d[d(g a Py/=g)/dt] i 



g a/3 dr 



(37) 



where L denotes the usual gravitational lagrangian density (eq.(10) in 
Although Tolman's formula was introduced as a measure of the total energy 
of the system, with no commitment to its localization, we shall define the 
mass within a sphere of radius r, inside S, as 



4tt f r V +A >/ 2 (T ° - Tl - 2T 2 2 ) dr 



1 

2 jo 



r 2 e (v+\)/2 



d 



dL 



dt \d[d(g a P-^=g) /dt] 



(38) 



is 



This (heuristic) extension of the global concept of energy to a local level 
suggested by the conspicuous role played by m T as the "effective gravitational 
mass" , which will be exhibited below. 

On the other hand, even though Tolman's definition is not without its prob- 



lems |22|, [23| , we shall see that my, as defined by (pq) , is a good measure of 
the active gravitational mass, at least for the system under consideration. 
After some simple but tedious calculations, it can be shown that (|3|) may 
be written as (see [23] for details). 



rax 



>+A)/2 



m(r, t) - Airr 6 T{ 



where the mass function m(r, t) is defined by |25|, |26 



m(r,t) = -rRl 32 



and the Riemann component for metric (|l]) is given by 

-^232 = 1 — e A 



(39) 



(40) 



(41) 
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Using field equations, ( ^0[ ) may be written in the most familiar form 

m{r,t) =4tt [ r 2 T°dr (42) 
Jo 

or alternativately |24| 

t) = ^ (T ° + Tl - T|) + r -C\, 2 (43) 



m(r 



where C| 32 denotes the corresponding component of the Weyl tensor. 
It is worth noticing that this is, formally, the same expression for in 
terms of m and Tl, that appears in the static (or quasi-static) case (eq.(25) 
in 13). 

Replacing Tl by (f|), and m by ( [40D and fl4"T|), one may also obtain 

mT = e ^- A )/Vy (44) 

This last equation brings out the physical meaning of as the active grav- 
itational mass. Indeed, it can be easily shown |27j that the gravitational 



acceleration (a) of a test particle, instantaneously at rest in a static gravi- 
tational field, as measured with standard rods and coordinate clock is given 
by 

(45) 



2 r 2 

A similar conclusion may be obtained by inspection of eq.(^) (valid only in 
the static or quasi-static case) P8|. In fact, the first term on the right side of 
this equation (the "gravitational force" term) is a product of the "passive" 
gravitational mass density (p + P r ) and a term proportional to m^/r 2 . 
We shall now consider another expression for my, which appears to be more 
suitable for the treatment of the problem under consideration. This latter 
expression will be evaluated immediately after the system departs from equi- 
librium. Therefore the physical meaning of mr as the active gravitational 
mass obtained for the static (and quasi-static) case, may be safely extrapo- 
lated to the non-static case within the time scale mentioned above. 
The required expression for the Tolman mass will be obtained as follows (see 
for details). Taking the r-derivative of (0) and using (|4*3]) and (|39D we 



obtain the following differential equation for tut 



11 



rm'rp — 3m T 



where W( s ) is given by 



>+A)/2 



e (X-u)/2 r 3 / 

+ —4— A+ 2 



4,rr 3 (r; - If) - 3W W 



(46) 



^,3g A 



1 i/'A' 



,'2 
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4 



4 



v 

T 



Equation 



can be formally integrated to obtain 



A' z/ N 
2r + 2r , 



(47) 



(m T )s — 



— r 



> / rE e (^+A)/2 



»7T 

r 



(T±-Tl) + ±J\itr\T«)>dr 



dr 



e (X-u)/2 / ^2 j^A 

- r 3 / AH } dr 

2r V 2 2 / 



(48) 



In the static (or quasi-static) case (A = A 2 = Az> = 0) the expression above 
is identical to eq.(32) in p0[ . 

We shall now proceed to evaluate (f4"8"[) immediately after perturbation. Using 
(|7|)-(|10"|) and (p7|)-(^0|), we see that up to first order we get immediately after 
perturbation (for both initial conditions) 



1 o 



■ P 
z>A 



inre 



Tl = 



(f+A)/2 



7^2 
J 2 



(p + P r ) uo + Qe 



Pi 



A/2 



Replacing (f49"D into ([||) and using ([33|) and (|36"D , we obtain finally 



(49) 



+ 4vrr 3 e {u+X)/2 \- (P r ~ P±) ~ \ f r^p'df 
Jr lr r 4 jo 



dr 



+ 4vrr 3 f E e A (p + P r ) w (1 - a) dr 



(50) 
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where the general expression for (?7Jt) s can be obtained from 
and (|16|) 



(m T ) s = m s + 



4?rr|g s (1 + 2cj s ) 
1-^ 



puj 



which, after perturbation reduces to 

(m T ) E = m E + 47rr|g s 



(O) 

(51) 

(52) 



6 Discussion 

Let us now consider a sphere of radius r within S. Immediately after pertur- 
bation the Tolman mass of this internal core is given by fl50|) . The relevance 



of the two terms in the first integral has already been discussed [23 and 
therefore they shall not be considered here. 

Instead, we shall focus on the last term in (|50|). If the system starts to 
collapse (Cj < 0) this last term tends to decrease the value of the Tolman 
mass, leading thereby to a weaker collapse. Inversely, if the system starts to 
expand (co > 0), the last term in fl5"0|) contributes positively to the Tolman 
mass of the core, leading to a weaker expansion. Thus, in both cases this 
last term tends to stabilize the system. This is so as long as a < 1. 
If a > 1 the inverse picture follows. In this case for initially collapsing 
(expanding) configurations the last term in fl5"0|) becomes positive (negative) 
leading to stronger collapse (expansion). 

In general the system becomes more and more unstable as a grows. 
It should be noticed that in the comments above we have assumed a to be 
constant througout the fluid distribution. This of course is a rather crude ap- 
proximation as it is evident from (p5[). Therefore, a wide variety of scenarios 
may be considered from different radial dependence of that parameter. 
Finally, observe that in the dissipationless case (a = 0) , an inflationary 
equation of state (p = — P r ) is equivalent to the critical point (a = 1) in the 
heat conducting situation. In both cases the stabilizer term in ( p0|) vanishes. 
In the heat conducting case a ^ 0, an inflationary equation of state leads to 
an effective inertial mass density equal to —kT/t, as follows from fl3"T| ) and 
(|33|). In this case, according to (|33|), ([48]) and (|49D the last integral in (|50|) 
should be replaced by 
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-4irr 3 / e A — udr (53) 

Jr T 

This last integral contributes negatively to Tolman mass in the case u > 0, 
yielding stronger expansions. 

In other words, in what concerns eq. (|50"D , an equation of state of the above 
mentioned form (p = — P r ) is equivalent (in the dissipative case) to a situation 
with a > 1. 

Of course one might ask if a real physical system may reach (or even go 
beyond) the critical point. The answer to this question seems to be affir- 
mative as suggested by the example provided in ||. Indeed it is shown in 
that reference that a mixture of matter and neutrinos with typical values of 
temperature and energy density, corresponding to the moment of birth of a 
neutron star in a supernova explosion may lead to values of a equal to or 
even greater than 1. 

However, it is not our purpose here to discuss about the plaussibility to reach 
the critical point but rather to bring out the physical meaning of a and the 
critical point. 

Finally it is worth noticing that evaluating the mass function from ( |3~9"D and 
using flSUP, we obtain similar conclusions about the relation between a and 
the mass function as those obtained for the Tolman mass. 
However, unlike the Tolman expression, the mass function can not be in- 
terpreted (for a part of the configuration) as the active gravitational mass 
and therefore the stability /instability criteria look less convincing when using 
m(r,t) instead of m^fr, t). 
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